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We discuss recently measured event-by-event fluctuations of transverse momentum and of mul- 
tiplicity in relativistic heavy-ion collisions. It is shown that the non-monotonic behavior of the 
Pt— fluctuations as a function of collision centrality can be fully explained by the observed non- 
monotonic multiplicity fluctuations. A possible mechanism responsible for the multiplicity fluctua- 
tions is also considered. 
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Event-by-event fluctuations of transverse momentum 
in heavy-ion collisions have been recently measured both 
at CERN SPS [1-3] and BNL RHIC [4-8], see also a brief 
review [9]. The data, which show a non-trivial behavior 
as a function of collision centrality, have been theoreti- 
cally discussed from very different points of view [10-27], 
including complete or partial equilibration [12,13,20,22], 
critical phenomena [14,27], string or cluster percolation 
[23,25], production of jets [11,26]. In spite of these ef- 
forts, a mechanism responsible for the fluctuations is 
far from being uniquely identified. Recently, the NA49 
Collaboration published very first data on multiplicity 
fluctuations as a function of collision centrality [28,29]. 
Unexpectedly, the ratio Var(N)/(N), where Var(N) is 
the variance and (N) is the average multiplicity of neg- 
ative particles, changes non-monotonically when number 
of wounded nucleons* grows. It is close to unity at fully 
peripheral (N w < 10) and completely central (N w > 250) 
collisions but it manifests a prominent peak at N w « 70, 
as shown in Fig. la. The measurement has been per- 
formed at the collision energy 158-A-GeV in the trans- 
verse momentum and pion rapidity intervals (0.005, 1.5) 
GeV and (4.0,5.5), respectively. The azimuthal accep- 
tance has been also limited, and about 20% of all pro- 
duced negative particles have been used in the analysis. 

The aim of this note is to show that the non-trivial be- 
havior of transverse momentum fluctuations can be ex- 
plained by the multiplicity fluctuations which enter the 
measures of pt— fluctuations. Specifically, we assume 
that in nucleus-nucleus collisions the event's transverse 
momentum is correlated to event's multiplicity exactly 
as in the proton-proton interactions [30] , and we express 
the <&— measure [31] of pt— fluctuations through the mul- 
tiplicity fluctuations. It is convenient for our discussion 



*A nucleon is called wounded if it interacts at least once 
in the course of nucleus-nucleus collision. The number of 
wounded nucleons N w approximately equals the number of 
participants, and we assume here that the equality holds. 



to use data on the transverse momentum and multiplic- 
ity fluctuations measured in the same experimental con- 
ditions. For this reason, we choose the data obtained by 
the NA49 collaboration which used the $— measure [31] 
to quantify the fluctuations of transverse momentum. 

Let us first introduce the measure. One defines the 
single-particle variable z = x — x with the overline de- 
noting averaging over a single particle inclusive distri- 
bution. Here, we identify x with the particle transverse 
momentum px- The event variable Z, which is a mul- 

tiparticle analog of z, is defined as Z =' YliLii^i ~ ^Oj 
where the summation runs over particles from a given 
event. By construction, (Z) = where (...) represents 
averaging over events. Finally, the measure is defined 
in the following way 
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It is evident that $ = 0, when no inter-particle correla- 
tions are present. Consequently, $ is 'deaf to the sta- 
tistical noise. The measure also possesses a less trivial 
property. Namely, <3 ) is independent of the distribution 
of number of particle sources if the sources are identi- 
cal and independent from each other [31,32]. Thus, the 
$— measure is 'blind' to the impact parameter variation 
as long as the 'physics' does not change with the collision 
centrality. In particular, the $ is independent of the im- 
pact parameter if the nucleus-nucleus collision is a simple 
superposition of nucleon- nucleon interactions. 

$(pt) measured in nucleus- nucleus collisions at SPS 
energy as a function of centrality [2] is shown in Fig. lb. 
The measurement has been performed in exactly the 
same experimental conditions as that of multiplicity fluc- 
tuations shown in Fig. la. As seen, both transverse mo- 
mentum fluctuations expressed in terms of <!> and multi- 
plicity fluctuations display a very similar centrality de- 
pendence, suggesting that they are related to each other. 

In the very first paper, where the <f>— measure was in- 
troduced [31], it was argued that the correlation between 
the event's multiplicity and transverse momentum is a 
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main source of the pt— fluctuations as quantified by $. 
For the case of p-p interactions, the problem was then 
studied in detail in [18]. Following this paper, we in- 
troduce the correlation (pt) vs. N through the mul- 
tiplicity dependent temperature or slope parameter of 
distribution. Specifically the single particle trans- 
verse momentum distribution in the events of multiplicity 
A is chosen in the form suggested by the thermal model 
i.e. 



P(n)(pt) ~ prexp 
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where m is the particle mass while TV is the multiplicity 
dependent temperature. In Ref. [18], XV was defined as 



TV = T + AT((A) - A) 



(2) 



with AT controlling the correlation strength. The param- 
eterization (2) was reasonable for proton-proton collisions 
where (A) is fixed, but it is not reasonable to study the 
centrality dependence in A-A collisions where (A) varies. 




FIG. 1. Multiplicity (a) and transverse momentum (b) fluc- 
tuations of negative particles as a function of number of 
wounded nucleons. The triangles correspond to p-p colli- 
sions, asterisks to C-C, squares to Si-Si, and circles to Pb-Pb. 
There are denoted statistical errors with vertical bars and to- 
tal errors including systematic uncertainties with horizontal 
dashes. The multiplicity data are taken from [28,29] while 
those on the transverse momentum from [2]. 



The correlation (px) vs. A at SPS energy, which is di- 
rectly observed [30,2] in p-p collisions, is most probably of 
simple kinematical origin. Namely, when the multiplicity 
of produced particles grows at fixed collision energy, there 
is less and less energy to be distributed among transverse 
degrees of freedom of produced particles. Consequently, 
the average event's pr decreases when TV grows. We ex- 
pect that the correlation (pr) vs. A is also present in A-A 
collisions at fixed centrality as the number of wounded 
nucleons controls the amount of energy to be used for 
particle production. However, we replace the parameter- 
ization (2) by 
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with ST = AT(N). The relation (3) correlates the slope 
parameter to the event's multiplicity A at fixed (A). 
The parameters T and ST are assumed to be indepen- 
dent of the centrality while the average multiplicity (A) 
depends (roughly linearly) on N w . As will be seen in our 
final formula (6), a small variation of T with the central- 
ity does not much matter. 

The inclusive transverse momentum distribution, 
which determines z 2 = p T — pr 2 , reads 
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where TV is the multiplicity distribution. The 
A— particle transverse momentum distribution in the 
events of multiplicity A is assumed to be the A— product 
of P(n){pt)- Therefore, all inter-particle correlations dif- 
ferent than (pr) vs. A are neglected here. Then, one 
easily finds 
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Assuming that the particles are massless and the corre- 
lation is weak, i.e. T ^> ST, the calculation of $ can be 
performed analytically. The result is [18] 



$(p T ) = V2 



{STf 



T(A) 5 
- (A 3 ) (A) 2 



((A 4 )(A) 2 -2(A 3 )(A 2 )(A) (4) 
(A 2 ) 3 + (A 2 ) 2 (A)) , 



where terms of the third and higher powers of ST have 
been neglected. As seen, the lowest non-vanishing con- 
tribution to <& is of the second order in ST. 

We intend to express <&(pt) through Var(N)/ (N) but 
&(pt), as given by Eq. (4), also depends on the third 
and fourth moment of multiplicity distribution. It would 
be in the spirit of our minimalist approach to use the 
multiplicity distribution which maximizes the Shannon's 
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information entropy S = Vn hiP/v [33] with (N) and 
Var(N) being fixed. The least biased method to obtain 
statistical distribution was prompted by Jaynes [34] . An 
application of the information theory to phenomenology 
of high-energy collisions is discussed in [35]. The multi- 
plicity distribution, which maximizes the entropy at fixed 
(A) and Var(N), is given by the formula 



V N = exp(a + bN + c A 2 ) , 



(5) 



where the parameters a, b and c are determined by the 
equations 

N N 

^(A - (N)) 2 V N = Var(N) . 

N 

Unfortunately there are no simple analytic expressions 
of a, b and c, and consequently the distribution (5) is 
very inconvenient to use. However, under the condition 
(A) > y / Var(N) > 1, which is usually satisfied in A- 
A collisions at fixed centrality, the distribution (5) can 
be replaced by the continuous Gauss distribution. Then, 
we get the required relations: ((A — (A)) 3 ) = and 
((A-(A)) 4 H3((V-(A)) 2 ) 2 . 
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FIG. 2. The average transverse momentum of negatively 
charged particles produced in p-p collisions as a function of 
event's negative particle multiplicity divided by the mean. 
The data are taken from [2] where the acceptance is precisely 
defined. The line corresponds to T = 137 MeV and 8T = 15.5 
MeV. 



Using these relations, the expression (4) gets the form 

r (STf Var(N) 
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Taking into account the already adopted assumption that 
(AO > y/Var(N) > 1, we finally find 



(§T) 2 Var(N) 



(6) 



When the Negative Binomial Distribution, instead of 
the Gauss, is used to describe the multiplicity distri- 
bution, one obtains the formula, which in the limit 
(A) > ^Var(A) > 1, coincides with Eq. (6). 

The values of the parameters T and ST for p-p col- 
lisions can be obtained from the NA49 data published 
in [2]. Following [18], we have computed the average pt 
at fixed N, using the distribution (1) with Tjv given by 
Eq. (3) . Comparing the results of our calculations with 
the experimental data [2], which are shown in Fig. 2, 
we have found T = 137 MeV and ST = 15.5 MeV. We 
note that T and ST are essentially independent from each 
other when the experimental data are fitted as ST deter- 
mines the slope of the curve shown in Fig. 2 while T 
controls its vertical position. For T = 137 MeV and 
ST = 15.5 MeV, the coefficient in the formula (6) equals 
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2.48 MeV 



(7) 



In Fig. 3 we compare the experimental values of <&(pt) 
with the predictions of the formula (6) with the numer- 
ical coefficient given by Eq. (7). As seen, the agreement 
is quite satisfactory 1 '. However, the analytic result (6) 
has been derived under several rather rough approxima- 
tions. Therefore, we have also performed a Monte Carlo 
simulation which is free of these approximations. For 
every nucleus-nucleus collision at a given centrality, we 
have first generated its multiplicity, using the Negative 
Binomial Distribution with the mean value and variance 
as in the experimental data [28,29]. Further, we have 
attributed the transverse momentum from the distribu- 
tion (1) with T = 137 MeV and ST = 15.5 MeV to each 
particle assuming, as in the experimental analysis [2], 
that all particles are pions. Having a sample of events 
for every centrality, the <&— measure has been computed. 
The statistical errors have been determined by means of 
the subsamplc method. The results of our simulation are 



Tn Ref. [18] the parameters T and 5T were estimated as 
167 and 8.2 MeV, using the data on p-p collisions at 205 GeV 
[30]. The data [2] were not available at that time. Then, 
V2(STY/T w 0.57, and the value of $(p T ) calculated by 
means of the formula (6) is dramatically underestimated. It 
was also concluded in [18] that the pr vs. N correlation pro- 
duces too small value of $(pt) to explain the experimental 
value. Now, this conclusion must be revoked. The discrep- 
ancy between the data [30] and [2] can be easily explained. 
The data [30] were collected at higher collision energy in the 
full phase-space while the NA49 measurement [2] was per- 
formed, as already mentioned, in the forward rapidity window 
(4, 5.5). 
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confronted with the experimental data in Fig. 3. As seen, 
there is a perfect agreement. 

Our calculations explicitly take into account only the 
Pt vs. N correlation. However, other correlations, in 
particular those due to quantum statistics, are not en- 
tirely neglected. Since we use the experimental value 
of Var(N)/(N), all correlations which contribute to this 
quantity, are effectively taken into account in our esti- 
mate of $(pt)- 

The multiplicity fluctuations strongly influence the 
Pt— fluctuations expressed in terms of <&, as the mea- 
sure $ depends on the particle multiplicity distribution, 
ft should be stressed that other fluctuations measures, 
such as F used by PHENIX Collaboration [4,5] or Act 
and (Tdyn by STAR [6,7], are also influenced by multiplic- 
ity fluctuations. Therefore, our main result (6) can be 
easily translated for F, Act or <7d yn - 
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FIG. 3. $(jjt) as a function of number of wounded nucle- 
ons. The open circles correspond to the NA49 data [2], the 
full circles show the results of our simulation while the trian- 
gles present the prediction of the analytical formula (6) with 
the numerical coefficient given by Eq. (7). 

A specific pattern of the pr— fluctuations have been ex- 
plained by the observed multiplicity fluctuations. Before 
closing our considerations we briefly consider a possible 
origin of the non-monotonic dependence of Var(N)/(N) 
on the collision centrality. For this purpose we first de- 
rive a well-known formula which relates particle number 
fluctuations to inter-particle correlations. The average 
multiplicity can be written as 



(AO 



/ d\p{v). 
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V is the volume of the interaction zone (fireball), where 
the particles are produced, and p(r) is the particle den- 
sity. The average multiplicity of produced particles is 
known to be roughly proportional to the number of 
wounded nucleons N w [36]. Since N w is in turn pro- 
portional to the volume V, we have (N) — p V with p 



being constant. The second moment of the multiplicity 
distribution can be written as 

(N(N-1))= f d 3 n [ d 3 r 2 p 2 (r 1 ,r 2 ) , 
Jv Jv 

where p 2 (ri,r 2 ) is the two-particle density. Defining the 
correlation function v(ri — r 2 ) through the equation 

P2(ri,r 2 ) = p(ri) p(r 2 )(l + ufa - r 2 )) , 

we get the desired formula 
Var(N) 



= l + p d 6 r v{y) 



(8) 



(N) j v 

which tells us that the multiplicity distribution is pois- 
sonian if particles are independent from each other (v(r) 
= 0). The pattern seen in Fig. la clearly shows that the 
particles are correlated at the stage of production. 

For further discussion we assume that the fireball is 
spherically symmetric and that its radius equals R ~ 

1 /3 

roNJ with ro^l fm. Then, the formula (8) reads 

= l + 47rp / drr 2 v(r) . (9) 
Jo 



Var(N) 



(N) 

It is not difficult to invent a correlation function v{r) 
which substituted in Eq. (9) reproduces the data shown 
in Fig. la. Various functions are discussed in [37]. Here 
we only describe qualitative features of v{r). The correla- 
tion function has to be positive at small distances (attrac- 
tive interaction) and negative at larger ones (repulsive 
interaction). The sign of the correlation changes at r k 
4 fm which corresponds to N w k, 70 when Var(N)/ (N) 
reaches its maximum. For r > (300) 1 / 3 w 7 fm the 
correlation function vanishes. A physical mechanism re- 
sponsible for such a correlation function is rather un- 
clear but some possibilities, which include: combination 
of strong and electromagnetic interaction, percolation, 
dipole-dipole interaction, and non-extensive thermody- 
namics, are discussed in [37]. 

We conclude our considerations as follows. A non- 
trivial behavior of transverse momentum fluctuations as 
a function of collision centrality can be fully explained 
by the centrality dependence of multiplicity fluctuations 
if the mean transverse momentum is correlated to the 
particle multiplicity in nucleus- nucleus collisions as in 
the proton-proton interactions. This correlation is most 
probably of simple kinematic origin. Our observation 
seem to exclude various exotic explanations of transverse 
momentum fluctuations. However, a mechanism respon- 
sible for multiplicity fluctuations still needs to be clari- 
fied. 
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